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Abstract 



In this letter we consider a prototype model which is described as an autonomous con- 
tinuous time delayed differential equation with just one variable. The chaos has been 
investigated with variable delay time and the synchronization phenomenon is examined 
both numerically and analytically using the Krasovskii-Lyapunov functions. We have ap- 
plied adaptive coupling law for synchronization, where the coupling equation also contains 
delay with modulated time. We also studied the effect of cryptography for this coupled 
system and the message extraction procedure is illustrated with the help of simulated 
results. 

1 Introduction: 

Synchronization between two dynamical systems [1] have stimulated a wide range of re- 
search activity [2] . The phenomena of synchronization in coupled systems have been es- 
pecially extensively studied in the context of laser dynamics, electronic circuits, chemical 
and biological systems [2]. Application of chaos synchronization can be found in secure 
communication, optimization of nonlinear system performance, modeling brain activity 
and pattern recognition phenomena [2]. The most important fields where chaos could 
have been practically applied is cryptography. The major concern in this field is that an 
encoded message is vulnerable for extraction by nonlinear dynamic forecasting, when it is 
hidden by the signal from the low-dimensional chaotic system. Then it becomes essential 
to develop high-dimensional chaotic systems having multiple number of positive Lyapunov 
exponents to implement a secure communication system. In this regard, one time-delayed 
system could be of a lot of attention: that is, 



For such systems with large constant time-delay r = tq [3] , some well-established nonlinear 
time-series analysis methods [4, 5, 6] ran into severe problems [7, 8]. 
In this letter we consider a simple one- variable prototype delay system [] containing cubic 
nonlinear ity. The system has a chaotic behavior for a small constant time delay which was 
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investigated in []. Here we consider some chaotic properties of this system with variable 
time delay. We also investigate the synchronization between two coupled systems with 
modulated time delay. The numerical results are supported by analytic calculation of the 
condition of chaos synchronization for both scalar and variable delay time with the help 
of Lyapunov functions. 

It is known that very simple time-delay systems are able to exhibit hyperchaos [9] . Conse- 
quently it has also been reported that time-delay systems provide alternative simple and 
efficient tools for chaos communication with low detectability []. Chaotic attractors of 
time-delay systems can have much higher dimension and many more positive Lyapunov 
exponents than the ordinary dynamical systems. 

In this letter, we also show that a message masked by chaos of a time-delay system with 
very high dimension and many positive Lyapunov exponents can be extracted successfully. 

2 Prototype delay model 

The delay differential equation describing the system are 



where 5 and e are positive system parameters, r is time delay. This system is used as 
a prototype model to observe self-oscillations in the shipbuilding industry. The chaotic 
behavior of system (2.1) for constant time delay are studied by A. Ucar []. Here we replace 
time delay parameter r as a function of time instead of constant delay as 



where tq is the zero frequency component, a is the amplitude and lo/2tt is the frequency of 
the modulation. The system are in chaotic state for the parameter value 5 = l,e = l,r = 
1.6. The behavior of the system can be change for variable delay time. For fixed value 
of ro = 1.6 and ujq = 0.8, we chose the amplitude parameter 'a' as bifurcation parameter. 
The system bifurcation diagram is depicted in figure (la), in which the system output x 
stereo-scopically observed from Poincare section and depicted versus the selected range of 
positive amplitude 'a'. Figure (la) shows the trajectory of the system has initially settled 
down at chaotic state for a < 0.275. For the range a e(0.275, 0.294) the system are in 
periodic states. With increase value of a, the system again reach in chaotic state. Further 
incresing a(> 0.3475) the system leads to unbounded solution namely unstable behavior. 
When the parameter 5 are change then the above bifurcation are changed. The parameter 
region between S and a are depicted in figure (lb). The system trajectory for a = 0.26 
and 5 = 1 are depicted in figure (lc). In the next section we consider the couple system 
and calculate analytically the sufficient condition for chaos synchronization. 

3 Coupled system and stability condition 

The couple system can be considered as 



x = 5x(t - t) - e[x(t - r)] 3 



(2.1) 



r(t) = tq + asin(ut) 



(2.2) 



x = g(x T , 5, e) = 5x(t — r) — e[x(t — r)] 3 
V = 9(Vt, 8, e) + k(x -y) = Sy(t - r) - e[y(t - r)] 3 + k(x - y) 



(3.1) 



where k is the coupling strength. 
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3.1 scalar delay time 

In this section, we study the sufficient condition for synchronization with the help of 
Krasovskii-Lyapunov theory. The desired synchronization manifold is expressed as x = y. 
Let A = x — y be the synchronization error. Then the dynamics of synchronization error 
is 

A = -kA + g'(x T ,S,e)A T (3.2) 

where g is the derivative of g with respect to time and A r = A(t — t). It is obvious 
that A = is the trivial solution of equation (3.2) for any value of time delay r. To study 
the stability of synchronization manifold x = y, we can use the Krasovskii-Lyapunov 
functional approach. According to K. PyragasfPRE, 58, 3067] the sufficient condition for 
stability of the trivial solution A = for the time delayed system 

A = -r{t)A + s(t)A T (3.3) 

is r(t) >| s(t) |. The Krasovskii-Lyapunov functional ( similar to Lyapunov function in 
the case of ODE) is 

V{t) = ^+fiJ A 2 {t + 6)d6 

where fj, > is an arbitrary positive parameter. The solution A = is stable if the 
derivative of the functional V(t) along the equation (3.3), V(t) = —r(t)A 2 + s(t)AA T + 
/uA 2 — fiA 2 is negative. V will be negative if 4(r — fi) > s 2 and r > fi > 0. The asymptotic 
stability condition for A = is given for 

fl = [^)A and r (t)>\s(t)\ (3.4) 

For a particular problem, two cases may arise, for the first case s is constant and r(t) is 
variable and the second case r is constant and s(t) is variable. The first case arises in 
[Prigas] and in this case fi is constant. But for the general cases the stability condition 
(i = i s not always true because [i is a parameter and s(t) is a variable. In our case, 
r(t) = k and s(t) = g'(x T ,5,e), where s(t) is a variable. So we can not apply the above 
result. But the above problem can be removed if we define /j as a function of time and 
the derivative of fi can be considered in the expression of V. 

Suppose fi = ((t) > 0, then the Krasovskii-Lyapunov functional can be taken as [PRE, 
75, 037203,2007] 

V (t) = ±A 2 (t)+at) jV(i+eR 

Then 

V(t) = -r{t)A 2 + s{t) AA T + C(t)(A 2 - A 2 ) + ((t) f° A 2 (t + £K 



If C(t) < then we have 

V(t)<-[r(t)-s 2 (t)/AC(t)-at)]A 2 
We obtain the stability condition as 

r(t)-s 2 (t)/A((t)-C(t)>0 

i.e. r(t) > h(s, C) where h(s, Q = s 2 (t)/4((t) + ((t) 

For any function of s(t), h(s,Q is a function of £(i) and has an absolute minimum for 
((t) = | and h m i n (s,C) =| s(t) \. Thus h(s,Q >| s(t) \ for any s and ( > 0. The stability 
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condition for synchronization is r(t) >\ s(t) |. For system (3.2) the stability condition for 
synchronization manifold x = y can be written as 



k >| g (x T , 5, e) | (3-5) 



3.2 Variable delay time 



When the time delay of the couple systems is modulated then the sufficient condition for 
synchronization are changed. For variable delay time one of the derivative terms in the 
expression of V(t) will be included. We define a positive definite Lyapunov functional of 
the form 



V{t)= l -tf(t)+at) f A 2 (H?R 

I J-r(t) 



Then 



V(t) = AA + C(t) [° A 2 (i + Ode + C(*)[A 2 - A 2 + A 2 r' 

J-r(t) 



if for ((t) < then V(t) < -A 2 F(X,((t)) where F(X,C(t)) = r(t) - ((t) - s(t)X - 
[C(t)r'-C(t)]X 2 withX = ^. 

In order to show that V(t) < for all A and A T i.e. for all X so it is sufficient 
to show that 

Fmin ^ 0. The absolute minimum of F occurs at X — 2C,{t)\i t 1 ) 
Fmin = r(t) — ((t) — s ' ^ , . Thus the sufficient condition for synchronization is 

Again ip(((t)) is a function of ((t). The minimum value of ip(((t)) occurs at 

With V'miri - M " 



\*(t)\ 



2V 



1-T 



Finally we get the sufficient condition for synchronization is 

Note that in the case of constant delay r = vanishes. So the above results are 

satisfied for constant delay. 

4 Numerical simulation 

In this section, we confirm that the numerical simulation fully can support the analytical 
results presented above. For constant coupling parameter and constant delay parameter 
in the couple system (3.1) we obtain r(t) = fc=constant and | s(t) |=| g (x T ,5,e) \. The 
stability condition (3.5) takes the form k >| g (x T ,5,e) |, where the maximum values are 

taken on driving system. The absolute maximum value is ^r^/J at x T = Thus the 

analytical condition for chaos synchronization is 

(4.1) 

For the parameter value e = l,r = 1 the synchronization region in the (k,S) plane is 
shown in figure (2a). At 5 = 1, for k = 0.2 and k = 1.8 the synchronization and non- 
synchronization states are shown in figures (2b) and (2c) respectively. 
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For time delay modulation the analytical condition for synchronization is 

I s(t) I . , 25 IT 1 , An . 

For the set of parameter value 5 = l,e = l,a = 0.2, uj = 0.8 the driving system are in 
chaotic state(fig lc), the analytical criterion for chaotic synchronization is k > 0.42. We 
take k = 1, the chaotic synchronization is shown in figure (3a). 

When the coupling strength of the couple system (3.1) are varied as 

x = Sx(t - t) - e[x(t - r)] 3 

y = 5y(t - r) - e[y(t - r)] 3 + k(t)(x - y) (4.3) 

k(t) = V (x - yf[k* - k(t - t 2 )} 

where rj = 1.0 is adaptive gain and k* = 1.9 is the minimum coupling strength when 
coupling are constant and T2 is the coupling delay time between systems x and y. r 
is given by equation (2.2) and we take coupling delay T2 as modulated delay time i.e. 
T2 = a — f3t 2 where a and (3 are constants. For a = 2 and (3 = 1, the variation of 
x(t) — y(t) and k(t) are shown in figures (3b) and (3c) respectively. 



5 Cryptography using synchronized system 

It is easy to see from Fig- (3a) that the error vector w = x — y tends to zero after some 
period of time t > to and we have x = y after to- We now wish to utilize these two sets 
of synchronized chaotic systems for the communication using the cryptographic encoding. 
The sender uses system () and the receiver uses system (). They choose the values of the 
variables x and y respectively as the secret keys after some time say t=10.0. To take the 
values of the secret keys as integer, they choose [lOOOx] and [lOOOy] respectively. Since 
they consider the coupled chaotic systems with time delay, to be more secured they both 
agree to change the values of r after every eight unit of messages. They both increases 
the values of r by 0.1 at this interval. The data from x(t) picked up by sender and the 
corresponding secret keys are shown in TABLE-1. 

Let us first describe the procedure how this is generated. The actual message that a 
sender wants to send to the receiver is called plaintext. The plaintext and the corre- 
sponding ciphertext can be divided into message units. P.G.Vaidya et al showed that each 
unit of message is a single alphabet and they used 26 letter equivalents 0, 1, ... 25. The 
corresponding formula for the ciphertext message is 

c = p + kmod (26) (2) 

and decrypted message can be obtained by 

p = c — kmod (26) (3) 

In the letter we have a generalization of the above procedure. Suppose our message 
contains both letters as well as numbers. In references [] the receiver recovered "GOOD- 
MORNING" in state of " GOOD MORNING" . An improved version of this method where 
the unit of message (it may be a word or a line or a paragraph) is not only alphabet 
but also numbers, decimals, space between two wards or sentences. The corresponding 
are merely 0,1,. . . 25 but more than that.We assign the numeric number to 9 by 0, 1, 
....9 respectively. The blank space(gap) between two word is assigned by 10. 11 is the 
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decimal (or full stop between two sentences or any punctuation marks). The 26 letters 
from A to Z are assigned from 12 to 37 respectively. The following table will give the clear 
picture of the complete assignment. 

Table-I 



Number assigned 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 




37 


Unit Message 





1 


2 


3 


4 


5 


6 


7 


8 


9 






A 


B 




Z 



The corresponding formula can be written as 



c« = Pi + ki mod(38) 
Pi = Oi — h mod(38) 

where ki are the secret keys to musk the message. Corresponding to every message unit we 
use one and only one message key, which are randomly generated. For a complete message 
the secret keys are a series of numbers {k±, &2, • • • k n }. Actually the key kj hide and secure 
the message unit pj. 

Let us consider the message which is a sentence that contains numbers also. 

WE HAVE OVER 2700 EMPLOYEES. 
The following table shows the message units and the corresponding keys. 

Table-II 



Keys 


h 


k 2 


^3 


&4 


k 5 


fc 6 


k 7 


k 8 


kg 


fcio 


kn 


ki2 


ki3 


ku 


ki5 


Unit Message 


w 


e 




h 


a 


V 


e 




o 


V 


e 


r 




2 


7 



Keys 


ki6 


kn 


kis 


hg 


k20 


k2i 


k22 


&23 


k2A 


k25 


k26 


k27 


&28 


Unit Message 










e 


m 


P 


1 


O 


y 


e 


e 


S 





Table-3 and Table-4 respectively shows the actual messages, corresponding plaintext and 
ciphertext and the message recovered from the receivers end. 
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Table-3 Ciphertext sent by the sender. 



Time 


Time 


x(t) 


Keys 


Plaintext 


Ciphertext 


Lag n 


t 




k 


P 


c = p+k mod (38) 


1.6 


200.0 


0.9012.. 


901 


W (34) 


61 




201.0 


0.9265.. 


926 


E (16) 


23 




202.0 


0.9296.. 


929 


-(10) 


27 




203.0 


0.9327.. 


932 


H (19) 


29 




204.0 


0.9515.. 


951 


A (12) 


13 




205.0 


0.9472.. 


947 


V (33) 


68 




206.0 


0.8412.. 


841 


E (16) 


21 




207.0 


0.8281.. 


828 


-(10) 


25 


1.8 


208.0 


1.4001.. 


1400 


(26) 


42 




209.0 


1.3874.. 


1387 


V (33) 


34 




210.0 


1.3574.. 


1357 


E (16) 


43 




211.0 


1.3353.. 


1335 


R (29) 


34 




212.0 


1.3227.. 


1322 


-(10) 


25 




213.0 


1.3090.. 


1309 


2(2) 


19 




214.0 


1.1356.. 


1135 


7(7) 


40 




215.0 


1.0755.. 


1075 


0(0) 


11 


2.00 


216.0 


0.47380.. 


473 


0(0) 


17 




217.0 


0.48044.. 


480 


-(10) 


22 




218.0 


0.48616.. 


486 


E (16) 


31 




219.0 


0.48616.. 


489 


M (24) 


57 




220.0 


0.50019.. 


500 


P (27) 


30 




221.0 


0.50342.. 


503 


L (23) 


32 




222.0 


0.51972.. 


519 


(26) 


51 




223.0 


0.52531.. 


525 


Y (36) 


67 


2.20 


224.0 


0.96425.. 


964 


E (16) 


23 




225.0 


0.95570.. 


955 


E (16) 


21 




226.0 


0.93889.. 


938 


S (30) 


43 




227.0 


0.91450.. 


914 


• (11) 


12 
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Table-4 Plaintext recovered by the receiver. 



Time 


Time 


y(t) 


Keys 


Ciphertext 


plaintext 


Lag n 


t 




k 


P 


p = c-k mod (38) 


1.6 


200.0 


0.9012.. 


901 


61 


34 (W) 




201.0 


0.9265.. 


926 


23 


16 (E) 




202.0 


0.9296.. 


929 


27 


10 (-) 




203.0 


0.9327.. 


932 


29 


19 (H) 




204.0 


0.9515.. 


951 


13 


12 (A) 




205.0 


0.9472.. 


947 


68 


33 (V) 




206.0 


0.8412.. 


841 


21 


16 (E) 




207.0 


0.8281.. 


828 


25 


10 (-) 


1.8 


208.0 


1.4001.. 


1400 


42 


26 (0) 




209.0 


1.3874.. 


1387 


34 


33 (V) 




210.0 


1.3574.. 


1357 


43 


16 (E) 




211.0 


1.3353.. 


1335 


34 


29 (R) 




212.0 


1.3227.. 


1322 


25 


10 (-) 




213.0 


1.3090.. 


1309 


19 


2(2) 




214.0 


1.1356.. 


1135 


40 


7(7) 




215.0 


1.0755.. 


1075 


11 


0(0) 


2.00 


216.0 


0.47380.. 


473 


17 


0(0) 




217.0 


0.48044.. 


480 


22 


10 (-) 




218.0 


0.48616.. 


486 


31 


16 (E) 




219.0 


0.48616.. 


489 


57 


24 (M) 




220.0 


0.50019.. 


500 


30 


27 (P) 




221.0 


0.50342.. 


503 


32 


23 (L) 




222.0 


0.51972.. 


519 


51 


26 (0) 




223.0 


0.52531.. 


525 


67 


36 (Y) 


2.20 


224.0 


0.96425.. 


964 


23 


16 (E) 




225.0 


0.95570.. 


955 


21 


16 (E) 




226.0 


0.93889.. 


938 


43 


30 (S) 




227.0 


0.91450.. 


914 


12 


11 (.) 



6 Conclusion 

In this communication, we have studied the chaotic behavior of prototype delayed model 
with variable time delay. The detailed parameter region and bifurcation diagram are ob- 
tained when the time delay of the system are modulated. The synchronization between 
two unidirectionally coupled chaotic systems with fixed and variable time delay are ana- 
lyzed, we have also obtained the sufficient condition for chaos synchronization manifold by 
K-L function and it was verified by numerical simulation. In the last section we considered 
the case of variable coupling in the modulated delay situation and have shown how it can 
be utilized in the process of communication using cryptography. 
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9 Caption of figures 

Figure (la) : Bifurcation diagram with respect to the amplitude parameter a when other 
parameters are S = 1, e = 1, r = 1.6, u = 0.8. 

Figure (lb): Parameter region between 5 and a where gray region is period one, light 
black region period two, deep black region are period three and white represent chaotic 
region. 

Figure (lc) : Phase space in the plane (x,x(t)) when a = 0.26. 

Figure (2a): Parameter region between coupling parameter k and system parameter 5. 
Shaded portion represent synchronization region. 

Figure (2b): Time variation of the relative error w = x — y, which indicates the system 
are not synchronized. 

Figure (2c): Time series of x — y which indicates the system are synchronized. 

Figure (3a): Time variation of the relative error w where the delay parameter r is a 

function of time. 

Figure (3b): Time evolution of w corresponding to the equation (4.3). 
Figure (3c): Time variation of the coupling parameter k. 
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